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Abstract 

The damping-induced self-recovery phenomenon refers to the fundamental 
property of underactuated mechanical systems: if an unactuated cyclic variable 
is under a viscous damping-like force and the system starts from rest, then the 
cyclic variable will always move back to its initial condition as the actuated 
variables come to stop. The regular momentum conservation phenomenon can 
be viewed as the limit of the damping-induced self-recovery phenomenon in the 
sense that the self-recovery phenomenon disappears as the damping goes to 
zero. This paper generalizes the past result on damping-induced self-recovery 
for the case of a single unactuated cyclic variable to the case of multiple unac- 
tuated cyclic variables. We characterize a class of external forces that induce 
new conserved quantities, which we call the damping-induced momenta. The 
damping-induced momenta yield first-order asymptotically stable dynamics for 
the unactuated cyclic variables under some conditions, thereby inducing the 
self-recovery phenomenon. It is also shown that the viscous damping-like forces 
impose bounds on the range of trajectories of the unactuated cyclic variables. 
Two examples are presented to demonstrate the analytical discoveries: the pla- 
nar pendulum with gimbal actuators and the three-link planar manipulator on 
a horizontal plane. 

Keywords: Mechanical system, cyclic variable, viscous damping, self-recovery 
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1 Introduction 



Mechanical systems of kinematically coupled structures, or multibody systems, have 
a variety of engineering applications including robotic manipulators, manufacturing 
machines with articulated components, bipedal walking and musculoskeletal system 
modeling [5]. An interesting aspect of such systems is how they behave when cer- 
tain degrees of freedom are left unactuated. In dealing with such under actuated 
mechanical systems, we often explore their fundamental properties. For instance, 
if the unactuated variables are cyclic (i.e. do not appear in the Lagrangian), then 
the momentum map associated with these variables will be conserved due to sym- 
metry |7| • One of the most celebrated examples of the symmetry and conservation 
of momentum is the falling cat problem [H [HJ [9]: a cat always lands on its feet 
after released upside down from rest, by deforming its body while keeping a zero 
angular momentum. There are variants of the falling cat such as Elroy's beanie 
[6] and springboard divers [2|. Another popular example is the system of a human 
sitting on a rotating stool holding a wheel pp. Sitting on the rotating stool, one 
spins the wheel by his hand while holding it horizontally. A reaction torque will 
be created and initiate the rotating motion of the stool in the opposite direction. 
As long as the wheel is rotating, the stool keeps rotating. After some time, if the 
person applies a braking force halting the wheel spin, then the stool will also stop. 
All of these phenomena follow from angular momentum conservation, which holds 
in the absence of external forces. If there is an external force, then the momentum 
is not conserved any longer in general and the consequent motion of the system will 
be different from that in the case of momentum conservation. 

We are here mainly interested in the case where the external forces are linear 
in the velocity, i.e. viscous damping-like forces. This was well studied in p] for the 
case of one-dimensional symmetry, i.e., S 1 or R 1 symmetry. According to [1], in the 
stool-wheel system with viscous damping friction on the rotation axis of the stool, 
the stool does not keep rotating but meets a bound in its motion even if the person on 
the stool keeps spinning the wheel. The larger the angular velocity of the spin of the 
wheel is, the more angle the stool rotates by, but in the end the stool meets a bound 
in its motion, not being able to keep rotating, as long as the angular velocity of the 
wheel is bounded. This bound on the motion of the stool is called the damping- 
induced bound. Another phenomenon, which is perhaps more remarkable, in this 
stool-wheel system with the friction is the following. After some time, if the person 
on the stool stops the spinning of the wheel, then the stool does not stop. Instead, it 
asymptotically converges back to its initial position, making the same number of net 
rotations in the opposite direction. This phenomenon is called damping-induced self- 
recovery. From this example, one can see that the viscous damping-like force plays 
a role of a restoring force, but it is different from the spring force. To the knowledge 
of the authors, an example of damping- induced self-recovery and boundedness was 
first reported by Andy Ruina at a conference |10| . where he shows a video of his 
experiment and provides an intuitive proof of the phenomena. The damping-induced 
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self-recovery phenomenon, without damping-induced boundedness, is also explained 
in [3] for the case where the damping coefficient is constant, but the proof therein 
is not complete. Both phenomena of damping-induced self-recovery and damping- 
induced boundedness were rigorously and completely proved for the first time in pQ , 
where the damping coefficient is allowed to be a function of the cyclic variable and 
may take some negative values. The results in [I] as well as \10\ |3] are for the case 
of one-dimensional symmetry. 

This paper generalizes the results in pQ to the case of higher-dimensional Abelian 
symmetry. We consider mechanical systems with several unactuated cyclic variables 
that are subject to external forces linear in the velocity. The systems are allowed to 
have some actuated variables that are subject to control forces. We find conditions 
on the linear forces under which the damping-induced phenomena of self-recovery 
and boundedness occur for the cyclic variables. The analytical approach taken in 
this paper for the proof of these two phenomena is different from that in pQ. The 
total order of M was implicitly made use of in pQ, but in this paper we employ a 
Lyapunov-like function to prove the damping-induced self-recovery and boundedness 
for the multiple cyclic variables. We first characterize a class of damping-like forces 
that induce new conserved quantities which we call in this paper the damping-added 
momenta, and then construct a Lyapunov-like function using the damping-added 
momenta to prove the existence of the phenomena of damping-induced self-recovery 
and boundedness. To illustrate the main results, we consider two nontrivial examples 
of multibody systems that possess multiple cyclic variables: the planar pendulum 
with gimbal actuators and the three-link planar manipulator on a horizontal plane. 

2 Main Results 

In this paper, the norm || • |j denotes the Euclidean norm for vectors and the corre- 
sponding induced norm for matrices. For an invertible matrix A = (atj), the (i,j)-th 
entry of its inverse matrix is denoted by a 13 . For a symmetric matrix A we denote 
its positive semi-definiteness by A y 0. For two symmetric matrices A and B of the 
same size, A y B means A — B >z 0. 

2.1 Equations of Motion and Damping- Added Momenta 

Let Q = Qi x Q2 be an n-dimensional configuration space, where Q\ = W and Q2 
is a smooth manifold of dimension n — r. If a unit circle S 1 appears as a factor of 
Qi, it is replaced by M, which does not impose any restrictions on the description 
of the dynamics. Let q = (x, y) € W x Q2 denote coordinates for Q, where 

x = (x a ) = (x\...,x r ), y = W) = (y r+1 ,...,y n ) 

and 

q=(q*) = (q\...,q r ;q r+ \...,q n ) = (x\...,x r ;y r+ \...,y n ). 
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For notational convenience, the following three groups of indices are used in this 
paper: 

l,...,rj r + l,...,n . 

a,/3,7,... a,b,c,... 



i,j,k... 

Consider a mechanical system with the Lagrangian 

L(qA) = -m a px a x p + m aa x a y a + -m ab y a y b - V(q) 

which is the kinetic minus potential energy of the system. Here we follow the 
Einstein summation convention. It is assumed that the mass matrix 



m 



m a p m ab 
m al3 m ab 



is symmetric and positive definite. 

We make the following four assumptions on the system: 



Al) The variables x a, s are cyclic [7], i.e. 

dL 







dx c 

for all a = 1, . . . , r. 
A2) Controls u a 's are given in the directions of y a, s. 

A3) Each cyclic variable x a is under a general damping-like force (i.e. linear in the 
velocity) described as — k a p(x)xP, where each coefficient k a p(x) is a continu- 
ously differentiable function of x. 



A4) The coefficients k a p satisfy 
and 



K/3 = 


k/3a 


dk a p 




dx"< 


dx@ 



(1) 

(2) 



By Al) - A3), the equations of motion of the system are written as 

d dL B 

S^ = '^' « = l»-.r (3) 
d dL dL 

— — — — = u a , a = r + 1, . . . , n. (4) 

dtdy a dy a v ' 

In individual coordinates, © and (HJ) can be expressed as 

m a piP + m ab y b + [ij, a}q l q ] = -k a piP (5) 

dV 

m a pxf + m ab y b + [ij, ajq'q 3 + — = u a , (6) 
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where [ij, k] denotes the Christoffel symbol 

r • • j.1 _ I ( dmik i dm i k _ dmi i 

for i, j, k = 1, . . . ,n. 

If there were no external forces, i.e. /c aj g = for all a, (3, then the momenta 
-g^k's would be the first integrals of the system ([3]), but due to the external forces 
F a = —kapx^'s, -^k's are not conserved any more. We will here show that there 
are new conserved quantities associated with (|3|) in place of gp«'s. 

Lemma 2.1. Assumption A4) entails the following: 

1. For each a, there is a function h a : W — > R such that 

2. For each h a , there is a function U : R r — > R such that 

Namely, the damping coefficient matrix (k a p{x)) is the second-order derivative 
matrix of U . 

Proof. 1. By ([2]) and Ponicare's lemma, there is a function h a that satisfies (J7J). 

2. By ([T]) and ((7J, we have dh a /dx^ = k a p = kp a = dhp/dx a . Hence, by 
Ponicare's lemma, there is a function U : W — > R such that (|SJ) holds. 

□ 

Using Lemma |2. 11 we now show that the dynamics (|3|) have r first integrals. 
Theorem 2.2. Let /i Q 's be functions that satisfy Q). Then, the r functions 

or equivalently 

mapx 13 + m Qa y a + /i Q (10) 
/or a = 1, . . . , r, are t/ie /irst integrals of the dynamics 

Proof. Differentiating ([9]) with respect to t along the trajectory of the system gives 
d f dL , \ d <9L 9/t r 



\dx a 

due to ([3]) and ([7]). This completes the proof. □ 

The vector-valued map (g, q) h-> (J^ + shall be called the damping-added 
momentum map. 
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2.2 Damping-Induced Self- Recovery and Damping-Induced Bound- 

edness 



The first integrals in Q depend on the initial condition (x(0), y(0), x(0), y(0)). Tak- 
ing an arbitrary initial condition and letting \i a be the initial value of the correspond- 
ing first integral, J^- + h a , we have 

m a pxP + m aa y a + h a (x) = fi a (11) 

for all t £ R. By Lemma 12.11 there is a function U on R r such that ([8]) holds. By 
adding —fi a x a to it, we can obtain a function Un : W — > R such that 

dU,.(x) , . . . . 

-J^J. = h a (x)-iM a . (12) 

The conservation equations (jlip for the damping-added momenta can be regarded as 
first-order differential equations for the cyclic variables x a, s. The damping-induced 
self-recovery phenomenon is a direct consequence of asymptotic stability of these 
first-order dynamics under some conditions on the function U^. Let us make the fol- 
lowing assumptions on which constitute sufficient conditions for the self-recovery 
phenomenon to occur: 

A5) The function Un has a unique critical point, denoted x e , and it is a minimum 
point of Un. 

A6) There is a number 6\ > such that 

inf U^(x) > U^(x e ). 

\\x—x e \\>5i 



inf ||d£^(x)||>0, 



A7) There is a number 82 > such that 

inf 

\\x — X e \\>S2 

where dU, = (gf, . . . , 

Assumption A6) guarantees that there does not exist any sequence {x^} in M r with 
limfc.j.oo ||xfc|| = 00 such that Wm.^^U ^(xu) = U^(x e ). Likewise, assumption A7) 
guarantees that there does not exist any sequence {x^} in W with lirm^oo \\xk\\ = 00 
such that lim^oo \\dU^(x k )\\ = 0. 

We now state one of the two main theorems of this paper. 

Theorem 2.3 (Damping-Induced Self- Recovery). Suppose that controls u a (t)'s are 
chosen such that lim^oo y(t) = and there are numbers c\ > 0, C2 > and C3 > 
such that 

cil 1 (m aP (y(t))) ± c 2 I (13) 
\\(m aa (y(t)))\\ <c 3 (14) 

for all t > 0. Then, lim^oo x(t) = x e and lim^oo x(t) = 0. In particular, if the 
initial condition is such that x(0) = and y(0) = 0, then lim^oo = x(0). 
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Proof. By (J7J), integration of ([3]) with respect to t yields (jlip for some constants 
// Q 's, which can be written as 

ir a = -m a ^-m Q %s a y a (15) 

by (fT2l) , By adding a constant, we may assume that U^(x e ) = 0. For each s > 0, 
define an open set W s in W' by 

W s = {x G M r | C/ M (x) < s}. 

Take any e > 0. We will find a number T > such that — a? e 1 1 < e for all t >T. 
By A5) and A6), there is a sufficiently small S > such that 

x e Ws ^ \\x — x e \\ < e. (16) 

Replacing 5 by a smaller positive number if necessary, we can find, by A5) - A7), 
an ei > such that 

x£W s => WdU^x)]] >ei. (17) 
Since the matrix (m a g) is symmetric and positive definite, (|13p implies 

-/ ^ (m°^(y(*))) ^ -I. (18) 

C2 Ci 

Choose a number < ^ < 1 such that 

i / i ^ 2 



which is always possible since the left hand side of (|19p is continuous in I and is 
positive at i = 0. Since linn^oo y(i) = and the matrix (m aa (y(t))) is bounded, 
there is a T\ > such that 

||(m aa (y(t))y a (t))|| < £ £l (20) 
for all t>Ti. Whenever x(t) <£ W s for some t > T 1} we have, by ([IS]), dHJ), {18]) 



7 



and U2QJ), 



dt 



dx c 
dx a 



-m 



m a,3 mp a y c 



a/3 (SUp 1 

-m H — — m aa y 

dx a 2 y 



< 



< 



< 




^ (jj^ + \m pb A + \m^m aa y a m pb y h 



+ ^\\{m aa y a ) II 2 
4ci 

4ci 



(21) 



where the right-hand side of (f2~TI) is negative by (fT9j) . Hence, U^(t) decreases at 
least linearly in time as long as x(t) ^ Ws- By definition of Ws, there must exist 
T > Ti such that x(t) € W s for all t > T. Thus, - x e \\ < e for all t > T by 

(fTo]h Therefore, lim^oo x(t) = x e . By taking the limit of both sides of (fTo]) . we 
obtain limt^ 00 x(t) = since dU^(x e ) = by A5). In particular, if x(Q) = and 
y(0) = 0, then x e = x(0) by A5) and CE]), so lim^oo x(t) = x(0). □ 



Remark 2.4. T/ie result in Theorem \2.3\ is global. In order to get a local result, 
one has only to assume A5) - A7) in a neighborhood of x e and to assume that the 
trajectory x(t) stays in the neighborhood. 

As was discovered in the case with a single cyclic variable [1] , the viscous damping 
force not only induces self-recovery but also imposes a bound to the range of the 
cyclic variable. Such a boundedness property also holds for multiple cyclic variables 
as stated in the following theorem. 

Theorem 2.5 (Damping-Induced Boundedness). Suppose that the function sat- 
isfies 

lim 11^(^)11= 00 (22) 



and 



lim U^x) 

\\x\\ — >oo 



CO. 



(23) 



If controls u a (t) 's are chosen such that y(t) is bounded and there exist c\ > 0, C2 > 
and C3 > such that ( f73|) and [1~4\ ) hold for all t > 0, then x(t) is bounded. 
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Proof. Since y{t) is bounded by assumption, there is a number C4 > such that 
||y(^)|| — c 4 f° r an t >0. By (|22p . there is a number ^ > such that 

lis " Sell > i => II^WII > 1 + 

ci 

Suppose that x(i) is not bounded. Then, there are numbers < t\ < t 2 such that 
\\x(t) - a; e || > i for all t € [ti,t 2 ] and E/^a^ti)) < U^xfa)) by $2$. Then, for all 

< -±\\dU,(x(t))\\ 2 + —\\dUMt))\\ 

C 2 c\ 

-l\\dU^x(t))\\ + ^j ||dC^(x(t))|| 

1 

< — • 

C2 



Hence, 



r*2 ^77 ^ 
< U„(x(t 2 )) - U„(x{h)) = I -£{t)dt < (t 2 - h) < 0, 



ti 



dt C 2 



which is a contradiction. Therefore, x(t) is bounded. □ 
2.3 Diagonal Damping Force: A Special Case 

Suppose that there are continuous functions ki,...,k r : M — > M. such that the 
damping coefficients k a p are given as 

ku(x) = k 22 (x) = k 2 (x 2 ); ■■■ ; k rr (x) = k r (x r ); (24) 

k a p{x) = for a ^ /3. 

Notice that we do not assume that fc a (x Q )'s take only non-negative values though 
we call them damping coefficients for convenience. 
A function h a satisfying (J7J) is given by 

h a (x a ) = / k a (s)ds, (25) 
J 

where there is no summation over the index a. Given [i = . . . ,fj, r ) E M r , the 
function defined by 



r ,. x a 

a=l Jo 



)ds (26) 



satisfies (PT2l). 
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Corollary 2.6 (Damping-Induced Self-Recovery). Suppose that the functions k a 's 
given in \2J$ and the functions h a 's defined in (|25p satisfy the following: 

(i) For each a, the equation h a (s) — fi a = has a unique root, which is denoted 
by x°. 

(ii) For each a, k a (x^) > 0. 

(Hi) For each a, there is an open interval I containing x® such that 

inf \h a (s) — u a \ > 0. 

sGR\/ 

Then, the function defined in (|26p satisfies A5) - A7) such that the conclusions 
in Theorem \2.S\ hold true. 

Proof. Take any a between 1 and r. Since h' a (x°) = k a (x®) > 0, the function 
h a (s) — fi a is increasing over an open interval containing x" . Thus, by condition 
(i), we have h a (s) — /i Q > for all s > x" and h a (s) — /i Q < for all s < x". It 
is now straightforward to show that the function defined in (|26j) satisfies A5) - 
A7). □ 

Corollary 2.7 (Damping-Induced Boundedness). Suppose that the functions h a 
defined in (1251) satisfy the three conidtions in Corollary \2.6\ and the following con- 
dition: 

(iv) For each a, lim^oo h a (s) = oo and lmis._j.__oo h a (s) = — oo. 
Then, the conclusion in Theorem \2.5\ holds true. 

Proof. It is easy to show that given in (126p satisfies (122p and (I23p . □ 

2.4 Choice of Control for Damping-Induced Self-Recovery 

In Theorems 12.31 and 12.51 we have assumed that an appropriate control law exists 
to satisfy some conditions on the trajectory y(t). We now constructively show that 
such a control law exists. 

Lemma 2.8 ([llj). The control law 

u a = fa(y,x,y) + g ab (y)T b , (27) 
where r h 's are new control variables and 

fa = ([ij, a] - m aa [ij, (3 }m ali ) q % q 3 - m aa m a ^ k^x 1 + — , 



dy a ' 



9ab = m ab - m aa m bl3 m al3 , 
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transforms the system in f5j) and to the following system: 

BJJ 

x a = -m a ^-m^m pb y b (28) 
V a = r a , (29) 
where Un is a function that satisfies (j 1 2|) . 

Proof. Solve ([5]) for xr , substitute it into ©, and apply the control in (J2TJ) to obtain 
(f29|) . Since dh a /dt = k a px^ , integration of ([3]) with respect to t yields ([TT]) where 
^ Q is the value of the first integral. Solving ([TT]) for x a gives ([28]) . □ 

Equation (|29p implies that we have full control over the motion of y a variables. 
Hence, it is always possible to find a control law to satisfy the assumptions in 
Theorems 12.31 and 12 . 5( refer to pQ for more details on the method of choosing controls 
u a 's. For example, suppose that y^(t) = (y^(t)) is a reference trajectory that y(t) 
must follow. Then, we can choose the following control law for r a : 

r a = m) - cUy a - m)) - c o(y a - *£(*)) (30) 

with cf > and Cg > such that the tracking error e a (t) := y a (t) — y°[(t) obeys the 
following exponentially stable dynamics: 

e a + c\e a + c a Q e a = 0. 

Thus, y a {t) and y a (t) converge exponentially to the reference trajectory y%(t) and 
y°[(t), respectively, for each a = r + 1, . . . , n. 

3 Examples 

In this section, we take two examples of mechanical systems with multiple cyclic vari- 
ables to demonstrate the phenomena of damping-induced self-recovery and bound- 
edness. 

3.1 Planar Pendulum with Gimbal Actuators 

Consider the planar pendulum system in Fig. [TJ The planar motion of the base block 
is unactuated and constrained to move freely along x and y directions only. The 
pendulum rod is assumed to be actuated by the gimbal-like mechanism which can 
apply torques along 6\ and Qi rotational axes. Choose coordinates as (q 1 , q 2 ; q 3 , q ) = 
(x, y; 9i,92), where x and y are unactuated but subject to viscous damping-like forces 
and 9\ and 62 are actuated by controls (i.e. n = 4 and r = 2). The mass matrix M 
is given by 

/ M a + Mb + m —mrc l c 2 mrs 1 s 2 \ 
_ M a + m mrc 2 

—mr c x c 2 mr 2 c 2 

\ mrs 1 s 2 mrc 2 mr 2 J 
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Figure 1: The planar pendulum with gimbal actuators. 

where Cj and Sj denote the cos#j and sin#j, respectively and the parameters in M 
are listed in Table [TJ Since x and y do not show up in M and also are not affected 
by the gravity, they are cyclic variables. 



Table 1: Parameters for planar pendulum 





Parameter 


Value 


Unit 


M a 


Slider mass 


2 


[kg] 


M b 


Gantry bar mass 


3 


[kg] 


m 


Pendulum ball mass 


3 


[kg] 


r 


Rod length 


0.5 


M 



To demonstrate the self-recovery phenomenon in this case, we apply the control 
law in Lemma 12.81 to control 6\ and 62 such that they move from [rad] to ^ [rad] 
and j [rad], respectively, within 3 seconds; see Fig. [2j For the unactuated cyclic 
joints, x and y, we simulate two different cases of damping matrix K = (k a fs)i< a ./3<2- 

f3 0\ T , /5 + 2cosx 4 
Kl = {o 3J' K2= { 4 4 + 2cosy 

Notice that K\ and K2 are the second-order derivative matrices of U\(x,y) = 
I (x 2 + y 2 ) and U2(x,y) = \x 2 + 2{x + y) 2 — 2 cos x — 2 cosy, respectively. One 
can easily verify that both U\ and U2 satisfy assumptions A5) - A7) and equations 
(1221) and flUD - 

The self-recovery phenomenon for both cases is verified in Fig. [3l From Fig. 



3(a) and Fig. [21 one can see that x(t) and y(t) converge to the origin, which is the 



12 




Actuated Joints (Gimbal Joints) 









K 2 




4 6 
Time (sec) 



10 



Figure 2: Time trajectories of the actuated variables, 9\{t) and 6*2 (i) . 



initial condition, after 6i(t) and 9 2 {t) settle down. The planar motion is described 
m Fig. |3(b)[ where we can clearly see the pendulum base automatically returns to 
its initial position for both cases. 



3.2 Three-Link Manipulator on a Horizontal Plane 

Consider the three-link open chain manipulator moving on a horizontal plane (i.e. 
no gravity effect) in Fig. HI One can easily see that the first joint angle 6\ is a 
cyclic variable, which is always the case for planar kinematic chains whose motion 
is constrained to a horizontal plane. For this particular system in Fig. HI the third 
joint angle #3 is also cyclic since the center of mass of the third link is located on 
its axis of rotation. Thus, we anticipate the self-recovery effect will occur in these 
two joint variables. To be compatible with the notation introduced in Section [21 we 
choose the configuration variables as q 1 = 0\, q 2 = 63 and q 5 = 9 2 . Then, the mass 
matrix is written as 

fa + 2f3c 2 h 5 + f3c 2 \ 
M=\ h h h , 
\5 + f3c 2 I 3 S J 

where c 2 denotes cos 6 2 , and parameters a, f3 and 5 are given by 

a = h + h + h + mir? + m 2 (l\ + r\) + m z {l\ + f 2 ), 

P = h(m 2 r 2 + m 3 £ 2 ), 

S = I 2 + h + m 2 rj + m 3 £l, 
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o 1 ^ 




-0.05 r 




5^ 1 1 1 1 1 * 

-0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0.05 



(b) x — y plot. 

Figure 3: Trajectories of unactuated cycle variables, x(t) and yit). 
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Figure 4: A three-link manipulator on a horizontal plane with 9\ and #3 unactuated. 

where individual parameter values are listed in Table [2 We use the following damp- 
ing coefficient matrix: 



which is the second-order derivative matrix of the function U{6\, #3) = | (2# 2 + #|). 



Table 2: Parameters for the three-link manipulator 





Parameter 


Value 


Unit 


h 


Link 1 length 


0.5 


H 




Link 2 length 


0.5 


H 


ri 


Location of link 1 center of mass 


0.1 


H 


r-2 


Location of link 2 center of mass 


0.1 


H 


h 


Link 1 moment of inertia 


2 


[kg • m 2 ] 


h 


Link 2 moment of inertia 


2 


[kg • m 2 ] 


h 


Link 3 moment of inertia 


2 


[kg • m 2 ] 


mi 


Link 1 mass 


10 


[kg] 


m 2 


Link 2 mass 


10 


[kg] 


m 3 


Link 3 mass 


10 


[kg] 



Figure [5] plots the time trajectories of the joint angles. The top plot shows the 
controlled motion of the second joint 62 while the second and the third present the 
time trajectories of the two unactuated joints Q\ and #3, respectively. A tracking 
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controller is designed according to ([30]) such that the actuated second joint #2 makes 
sixteen full revolutions (i.e. 32tt [rad]) during the first 8 seconds of time. The initial 
values for the joint angles are set to 6*2(0) = § , and #i(0) = #3(0) = 0. Figured] 
clearly shows that both of the unactuated cyclic variables 9\ and #3 converge back 
to their initial positions right after O2 is regulated at 32tt [rad]. Note that 9\ returns 
back to its original position after having made two full turns, demonstrating that the 
self-recovery effect is global. Besides the self-recovery phenomenon, the trajectories 
of 9\ and #3 in Fig. [5] also demonstrate the damping-induced boundedness derived 
in Theorem 12.51 More specifically, we can observe that 6\ and 62 oscillate around 
— 4-7T and — 3ir, respectively while 62 is still increasing. 
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Figure 5: Time trajectories of joint angles, 9%(t) (actuated), 9\ (t) (unactuated), and 
#3 (t) (unactuated). 



4 Conclusion and Future Work 

The phenomena of damping-induced self-recovery and damping-induced bounded- 
ness have been generalized to mechanical systems with several unactuated cyclic 
variables. The major contribution of this paper comes from characterizing the 
damping coefficient matrix as the second-order derivative matrix of a function and 
identifying a class of such functions that guarantee the self-recovery and bounded- 
ness of all unactuated cyclic variables. Regular momentum conservation is the limit 
of the damping-induced self-recovery as the damping disappears, in the sense that 
the recovery phenomenon vanishes in the limit. Non-trivial examples of mechanical 
systems with multiple cyclic variables are provided to demonstrate the theoretical 
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discoveries. 

It will be interesting to study the possibility of occurrence of damping-induced 
self-recovery and boundedness for the case of non-Abelian symmetry in the unactu- 
ated variables. We are currently examining the system of a spacecraft with internal 
rotors, where the spacecraft experiences viscous damping friction as it rotates. 
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